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Abstract 

We analyze the formation of the low-energy structure (LES) in above-threshold ionization spectra 

n r 

first observed by Quan et al. [1| and Blaga et al. [2J using both quasi-classical and quantum 
approaches. We show this structure to be largely classical in origin resulting from a two-dimensional 
focusing in the energy- angular momentum plane of the strong-field dynamics in the presence of 
the atomic potential. The peak at low energy is strongly correlated with high angular momenta of 
the photoelectron. Quantum simulations confirm this scenario. Resulting parameter dependences 



agree with experimental findings [l|, |2j and, in part, with other simulations 
PACS numbers: 32.80.Rm,32.80.Fb 
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One of the important fundamental processes of strong-field laser-matter interaction is 



Ml- 



above-threshold ionization (ATI), a process generally considered to be well understood p- 
The workhorse for the description of laser-matter interaction, the strong-field approximation 
(SEA), is expected to be valid in the regime of small Keldysh parameters 7 = (I p /2U P ) 1 ^ 2 , 
where I p is the ionization potential (binding energy) of the electron to be ionized and 
Up = Fq/Alu 2 the ponderomotive energy associated with the free quiver motion in the laser 
field with amplitude a = F /u 2 . Small 7 < 1 corresponds to large field strength F and 
low frequency u. In this so-called tunne ing regime the SEA (and its semiclassical incarna- 
tion, the "simple-man's model" (SMM) 9j) should work well as the influence of the atomic 
Coulomb potential is reduced to a weak perturbation. It thus came as a major surprise 
when for short mid-infrared (mid-IR) laser pulses A > 1600 nm and intensities of I ~ 10 14 
W/cm 2 (7 < J) an unexpected peak-like low-energy structure (LES) was found in ATI 
spectra p], [2I, [h| in contradiction to the SFA featuring a smooth photoelectron spectrum. 
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Therefore, the LES was termed an "ionization surprise" 

First experimental studies revealed several characteristic features of the LES providing 
important clues as to its origin. The LES was found to be universally present irrespective 
of the ionized atomic or molecular target in the energy range between 1 to 10 eV, its peak 
position and width primarily depend on the laser parameters and are approximately pro- 
portional to the ponderomotive potential Up oc 7~ 2 , and the LES is absent for circularly 
polarized laser fields. The latter finding seems to indicate that re-scattering, more precisely, 
revisiting of the ionic core by the electron plays a significant role. Consequently, longitudi- 
nal (i.e. along the laser-polarization axis) and transverse Coulomb scattering were invoked 
in theoretical explanations of the LES which was found to be present in the full numerical 
solution of the time-dependent Schrodinger equation (TDSE) 2j. Meanwhile, several mod- 
els have been proposed invoking different and, in part, contradictory notions. While the 
coupling between the longitudinal and transverse degrees of freedom by Coulomb scattering 
was emphasized in 3|, |4| , the appearance of the LES was demonstrated in a strictly one- 
dimensional model in the absence of a transverse degree of freedom { 
LES feature was also found for screened (short-ranged) model potentials 
as to the role of long-range Coulomb scattering and focusing. 

In the present paper, we relate the LES to a classical two-dimensional focusing in phase 
space which gives rise to strong correlation between the energy E and the angular momentum 



Furthermore, an 
2| raising questions 



L of the LES electrons since for small 7 or large U p (compared to the discrete level spacing) 
a close classical-quantum correspondence is to be expected. We perform quasi-classical 
simulations using the classical trajectory Monte-Carlo (CTMC) method including tunneling 
Il2l which is validated by full TDSE simulations confirming the E—L correlation. We deduce 

nn 

laser-parameter dependences of the LES consistent with experiment 12[. Moreover, we 
identify a pronounced carrier-envelope phase (CEP) dependence of the LES which may open 
another route towards experimentally monitoring the CEP of mid-infrared pulses. Atomic 
units are used throughout this paper unless otherwise stated. 

Our CTMC simulation employs a standard adaptive step-size Runge-Kutta propagator 



with initial conditions chosen following 13J: the starting coordinate on the polarization axis 
(xq = yo = 0) is given by the tunnel exit in the combined Coulomb and laser fields, 

_ y Fp + y/(i p /F y - 4/F0 ^ 1 
z °~ 2 ~ F ' [ ) 

Both longitudinal (pn) and transverse (p±) momentum distributions at the tunnel exit are 



Gaussian distributed centered at p\\ t ± = 13J. While the initial longitudinal momentum 



distribution turns out to be unimportant for the LES (one could set pii = without significant 
change of the results), the width of the transverse momentum distribution 



P± = J-i= ( 2 ) 




is key to sampling the relevant phase-space region within which focusing of the LES occurs. 
Typically, 10 6 initial conditions from the distribution are sampled and are propagated in 
the combined laser and atomic Coulomb fields. At the end of the pulse the asymptotic 
momentum distribution is determined analytically by propagating along Kepler orbits [ijj]. 
For the latter we also allow exponential screening to probe for finite-range effects. 

The classical phase-space analysis is checked against full 3D quantum-dynamics simula- 
tions. The TDSE is solved by discretizing the coordinate space in a pseudospectral grid 

and prop agating the wave function by the split-operation method in the energy representa- 
I I 

tion [15|. When the time-dependent wave function in space reaches the outer region where 
the atomic potential becomes negligible compared to the kinetic energy, we project the 



outer region wave function on Volkov states to obtain the momentum distribution 16]. This 
method allows for long-time propagation of the wavepacket without encountering unphysical 
reflections at the boundary. 
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FIG. 1. (Color online) Distribution of ionized electrons in the energy (E) and angular momentum 
(L) plane observed within an angular cone of 6 C = ±10° around the polarization axis. The island 
at large values of L is the origin of the LES as demonstrated by the projection onto the energy 
axis. Inset: electric field of laser pulse (2200 nm, 10 14 W/cm 2 , 8 cycles) 

As one typical reference laser pulse we choose a cosine-like pulse shape (carrier-envelope 
phase 4>cep — 0) with sin 2 -envelope for the electric field, a wavelength of A = 2200 nm, a 
peak intensity of I = 10 14 W/cm 2 , and a duration of 8 cycles (total duration ~ 60 fs, full 
width at half maximum (FWHM) of intensity ~ 20 fs). Variation of the laser parameters 
will be discussed below. For this pulse shape tunneling ionization is strongly concentrated 
near the local field maxima of three adjacent half-cycle pulses at the center of the pulse. 
The final distribution in the energy (E) - angular momentum (L) plane (Fig. [1]) of electrons 
emitted into a double cone oriented along the polarization axis (0° and 180°) with opening 
angle 8 C = ±10° features an island with three distinct peaks at high angular momenta. Upon 
projection onto the energy axis, this high-L-low-i? island can be identified as the source of 
the LES with an energy of E ps 3 eV. This island is well-separated from the "background" 
of low-energy low-angular momentum electrons. The sharp parabolic boundary of the latter 
is given by 

E = 2a 2 sin 2 6 C • (3) 
The positions (Ei, Li) of the peaks within the island are determined by the variation of 
the subsequent field maxima within the pulse envelope. Note that they are not related to 
the peak sequence due to multiple reversals in a monochromatic (constant amplitude) field 
identified in the ID analysis [5]. A full quantum simulation for identical laser parameters 
(Fig. [2]) confirms the appearance of the high-L island. Note that because of the Heisen- 
berg uncertainty relation (ALA6) > 1 a straight-forward comparison between classical and 
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FIG. 2. (Color online) Comparison between classical (a) and quantum (b) E — L distributions for 
angle-integrated emission (same laser parameters as in Fig. [I]) . 

quantum E — L distributions is possible only for the angle-integrated emission. 

It is now instructive to explore within classical dynamics to the origin of the high-L-low- 
E island. Two prototypical trajectories reaching this island (Fig. [3]) visualize the strong-field 
dynamics perturbed by the atomic field (in the present example: a pure Coulomb field). Both 
trajectories are launched at the tunnel exit (z pa —10 a.u.) near the same field maximum 
with slightly different transverse momenta. The subsequent quiver motion is perturbed by 
the combined Coulomb and laser fields upon return to the vicinity of the ionic core. The 
trajectories eventually lock on to Kepler hyperbolae [ijj of different orientation but similarly 
large angular momenta. The rotation of the major axis of the Kepler hyperbola leading to 
emission into the lower hemisphere (dashed blue in Fig. [3]) is consistent with the semiclassical 
trajectories identified in (J). We emphasize that such an axis rotation is one pathway but 
not necessarily a prerequisite for reaching the LES as illustrated by the trajectory (solid red 
line in Fig. [3]) emitted into the same hemisphere as originally launched. The Coulomb field 
provides a "kick" when the slowing down of the quiver motion near its turning point occurs 
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FIG. 3. (Color online) Two typical trajectories contributing to the LES launched during the same 
half-cycle and eventually locked onto Kepler hyperbolae with the nucleus in its focus (thin lines). 
When the transverse momentum transfer near the turning point exceeds the initial transverse 
momentum the orientation of the Kepler hyperbola flips. 

in the vicinity of the ionic core. Such a scenario of the strong-field quiver motion only locally 
perturbed near the "inner" turning points requires a large quiver amplitude a compared to 
the typical distance from the core, taken to be of the order of the tunnel exit z , 

- - T7?V 2 * 4 » 1 ( 4 ) 
z I p {2ttc) 2 Y 

and readily explains why a pronounced LES develops only for mid-infrared pulses. The 
control parameter (Eq. H]) a/zQ = 2/7 2 yields the experimentally observed scaling of the 
upper border of the LES, Eh ~ 7 -1 ' 8 , and accounts for the overall variation of the LES with 
7 (but not with / and oo separately) 

The formation of the high-L island for the CTMC ensemble is the signature of a two- 
dimensional focusing in phase space. Key quantity is the classical transition probability 
\v\ for mapping the initial conditions, the time t (or laser phase (p = out), and the initial 
transverse momentum p± with which the trajectory is launched onto the final-state variables 
E and L given by the Jacobian of the dynamical transformation 

\d((f,p± 



P(E,L) 



(5) 



d(E,L) 

In Eq. [5] the range of p± is given by the ADK transverse momentum distribution and ip is 
the phase angle of the laser field relative to a given field extremum F j close to which the 
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FIG. 4. (Color online) Visualization of the Jacobian of the mapping (<p,p±) — > (E,L) for initial 
conditions near the absolute field maximum of a cosine-pulse at t = a.u. Shown is the dynamical 
mapping of the plaquettes (A(fj, Ap±j) j = 1...6 onto the E — L plane. Contours indicate 
the initial ADK ionization probability (top) and P(E,L) (bottom). To illustrate 2-dimensional 
mapping the longitudinal momentum is set to pu = a.u. and the position of the tunnel exit to 
zq = 7 a.u. 



launch of the trajectory occurs. The phase space focusing can be directly visualized by the 
distortion (stretching and compression) of the plaquettes (Aipj, Ap±j) of initial combinations 
when mapped onto the E — L plane (Fig. 0J. In close analogy to pulse-induced focusing of 
Rydberg wavepackets [3], focusing in the E — L plane proceeds via the confluence of two 
extrema in the deflection function to a point (line) of inflection. Traces of this 2D focusing 
are present also in the projection onto the energy axis. This explains the appearance of the 
LES feature and of focusing also in reduced- dimensional models j^]. The position of the 
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resulting peak (Ei, Li) can be estimated from 
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P±,t t~ ~ Px,tai+i ; ( 7 ) 



where Eq. [6] follows from the impulsive momentum transfer during the longitudinal velocity 
reversal at the turning point close to the nucleus in the (first) subsequent half-cycle while 
iCq. [7] follows from the locking of the quiver motion onto the asymptotic Kepler hyperbola 



14l | . Note that the final coordinates (Ei,Li) for electrons emitted during the ionization 
burst near the z-th field maximum are primarily determined by the field strength of the 
subsequent field maximum -Fcm+i. This observation explains the temporal order of the peaks 
within the island in E and L due to three subsequent ionization events (see Fig. [[]). A 
further consequence is that for longer pulses with a slowly varying envelope the structure of 
individual peaks within the island merge into a ridge whose upper border in both E and L 
is determined by the absolute field maximum at the center of the pulse. 

As the formation of the high-L island results from the interplay between the strong-field 
dynamics and the perturbation by the atomic force field near the inner turning point, the 
dependence on the range of the atomic potential is of conceptual interest. We therefore 



employ an exponentially screened potential 

V d {r) = -^e-n\ (S) 



r 

where d is the screening length of the potential and Z(d) is adjusted such that the ionization 
potential remains constant at its hydrogenic value (I p = 0.5 a.u.). The Coulomb limit 
corresponds to d = oo and Z(oo) = 1. With decreasing d the high-L structure moves for 
d = 10 towards lower L ~ 10 and disappears entirely for a truly short-ranged potential 
with d ~ 1 (Fig. |5]). Obviously, the formation of the LES-generating island requires atomic 
potentials whose range d extends to at least to a distance of the order of the tunnel exit 
(^o ~ d < a ] z o ~ 10 a.u. in the present case) such that the atomic field can impart a 
longitudinal momentum "kick" to accelerate and a transversal "kick" to deflect the outgoing 
trajectory at the inner turning point of the subsequent half-cycle. In agreement with [2j we 
find a pure Coulomb field (d = oo) is thus not necessarily a prerequisite for the occurrence 
of the LES, the shape and position of the high-L-low-E' island forming the LES, however, 
is very sensitive to the interactions at intermediate to large distances. 
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FIG. 5. (Color online) Variation of the angular integral E — L distribution with screening length d: 
top: d = oo, Z = 1 (Coulomb), center: d = 10, Z = 1.1, bottom: d = 0.9, Z = 2. Laser parameters 
as in Fig. [H observation angle 9 C = ±10°. The resulting electron spectra are shown as projection 
onto the energy axis. 

The focusing effect, identified here as underlying the LES, opens up novel opportunities 
to monitor the carrier-envelope phase of ultrashort few-cycle mid-infrared pulses. As tra- 
jectories asymptotically propagate in the direction opposite of that of the tunnel exit, we 
expect different spectra in forward (defined by the direction of the maximum laser field) and 
backward directions (Fig. [6]). For a short cosine-pulse (here 8 cycles; 4>cep = 0) the LES 
in both position and shape differs significantly from the sine-shaped pulse (4>cep = 7r/2, 
bottom panel). For the cosine-pulse, the emission at the absolute field maximum leads to a 
single peak near E = 2.5 eV and 6 ~ 180° observation angle while the two adjacent local 
maxima give a double peak structure at E = 2 and E = 3 eV at 9 m 0°. Such differences 
between forward and backward directions can also be observed for other few-cycle pulses. 
If the carrier-envelope phase 4>cep of the laser field is shifted by n/2 (8 cycles, sine-pulse, 
bottom panel) both forward and backward LES consist of only one peak. Due to the energy 
shift of the LES per half-cycle, these LES are observed at different energies (forward ~ 2.6 
eV, backward ~ 3 eV). This promises to give access to information on the CEP-phase 4>cep- 
To our knowledge, this is the first identification of signatures of carrier-envelope phase effects 
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FIG. 6. (Color online) Angle-resolved energy spectra of photoelectrons emitted in the interaction 
of hydrogen with an 8 cycle laser pulse with A = 2200 nm and I = 10 14 W/cm 2 . Top: cosine-pulse 
((frcEP = 0), bottom: sine-pulse (<Pcep = 7r /2) 

in the low-energy part of the photoelectron spectrum. Previous studies were focused on the 
high-energy tail of the spectrum and shorter wavelengths (A ~ 800 nm, |19|). 

In conclusion, we have shown that the recently discovered low-energy structure in strong- 
field ionization results from a two-dimensional focusing effect of the classical phase-space 
distribution at high angular momenta and low energy. Full 3D quantum simulations confirm 
this scenario. This structure appears for small Keldysh parameters 7 or, equivalently, large 
ratios of quiver amplitude a to distances of the tunnel exit z , a/z ^> 1. In this regime, the 
strong-field quiver dynamics is locally perturbed by the atomic force field near the "inner" 
turning point during the half-cycles subsequent to the one causing tunneling ionization. 
While the position and shape of the peak in the E — L plane is sensitively influenced by 
the long-range Coulomb interaction, the appearance of an enhancement occurs also for a 
screened atomic force field provided its range is close to or exceeds zq. We show the the LES 
for few-cycle pulses is sensitive to the carrier-envelope phase which may open the opportunity 
to measure the CEP in low-energy spectra. 
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